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Abstract 

The paper deals with a spatially homogeneous and isotropic FRW space-time filled with perfect 
fluid and dark energy components. The two sources are assumed to interact minimally, and therefore 
their energy momentum tensors are conserved separately. A special law of variation for the Hubble 
parameter proposed by Berman (1983) has been utilized to solve the field equations. The Berman's 
law yields two explicit forms of the scale factor governing the FRW space-time and constant values of 
deceleration parameter. The role of dark energy with variable equation of state parameter has been 
studied in detail in the evolution of FRW universe. It has been found that dark energy dominates 
' the universe at the present epoch, which is consistent with the observations. The physical behavior 

of the universe is discussed in detail. 
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1 Introduction 

Recent observations of type la supernovae (SN la) [1-5], galaxy redshift surveys [BJ, cosmic microwave 
background radiation (CMBR) data |8] and large scale structure j9] strongly suggest that the ob- 
servable universe is undergoing an accelerated expansion. Observations also suggest that there had 
been a transition of the universe from the earlier deceleration phase to the recent acceleration phase 
|10j . The cause of this sudden transition and the source of the accelerated expansion is still unknown. 



Measurements of CMBR anisotropics, most recently by the WMAP satellite, indicate that the universe 
is very close to flat. For a flat universe, its energy density must be equal to a certain critical density, 
which demands a huge contribution from some unknown energy stuff. Thus, the observational effects 
like the cosmic acceleration, sudden transition, flatness of universe and many more, need explanation. 
It is generally believed that some sort of 'dark energy' (DE) is pervading the whole universe. It is a 
hypothetical form of energy that permeates all of space and tends to increase the rate of expansion 
of the universe [TT] . The most recent WMAP observations indicate that DE accounts for 72% of the 
total mass energy of the universe [12] . However, the nature of DE is still a mystery. 

Many cosmologists believe that the simplest candidate for the DE is the cosmological constant 
(A) or vacuum energy since it fits the observational data well. During the cosmological evolution, 
the A-term has the constant energy density and pressure p( de ^ = —p( de \ where the superscript (de) 
stands for DE. However, one has the reason to dislike the cosmological constant since it always suffers 
from the theoretical problems such as the "fine-tuning" and "cosmic coincidence" puzzles [13]. That 
is why, the different forms of dynamically changing DE with an effective equation of state (EoS), 

u (de) = p (de) 

I p( de ) < —1/3, have been proposed in the literature. Other possible forms of DE include 
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quintessence (u/ de ) > — 1)[H], phantom (u/ de ) < — 1) [15] etc. While the possibility « — 1 is 

ruled out by current cosmological data from SN la (Supernovae Legacy Survey, Gold sample of Hubble 
Space Telescope) [U [16], CMBR (WMAP, BOOMERANG) pU EE] and large scale structure (Sloan 
Digital Sky Survey) [19J data, the dynamically evolving DE crossing the phantom divide line (PDL) 
(a/ rfe ) = —1) is mildly favored. SN la data collaborated with CMBR anisotropy and galaxy clustring 
statistics suggest that -1.33 < u/ de) < -0.79 (see, Tegmark et al. [20]). 

In 1983, Berman |21] proposed a special law of variation of Hubble parameter in FRW space-time, 
which yields a constant value of deceleration parameter (DP). Such a law of variation for Hubble's 
parameter is not inconsistent with the observations and is also approximately valid for slowly time- 
varying DP models. The law provides explicit forms of scale factors governing the FRW universe and 
facilitates to describe accelerating as well as decelerating modes of evolution of the universe. Models 
with constant DP have been extensively studied in the literature in different contexts (see, Kumar 
and Singh [22] and references therein). Most of the models with constant DP have been studied by 
considering perfect fluid or ordinary matter in the universe. But the ordinary matter is not enough to 
describe the dynamics of an accelerating universe as mentioned earlier. This motivates the researchers 
to consider the models of the universe filled with some exotic type of matter such as the DE along with 
the usual perfect fluid. Recently, some dark energy models with constant DP have been investigated 
by Akarsu and Kilinc [231 [Ml EH] , Yadav [26j and Yadav and Yadav [27] . 

In this paper, we have considered minimally interacting perfect fluid and DE energy components 
with constant DP within the framework of a FRW space-time in general relativity. The paper is 
organized as follows. In Sect. 2, the model and field equations have been presented. Sect. 3 deals 
with the exact solutions of the field equations and physical behavior of the model. Finally, concluding 
remarks have been given in Sect. 4. 



2 Model and field equations 

In standard spherical coordinates (x l ) = (t,r,8,(ft), a spatially homogeneous and isotropic FRW line 
element has the form (in units c = 1) 



ds 2 = -dt 2 + a 2 (t) 



(h2 +r 2 (^ 2 + sin 2 



1 — kr 2 



(1) 



where a(t) is the cosmic scale factor, which describes how the distances (scales) change in an expanding 
or contracting universe, and is related to the redshift of the 3-space; k is the curvature parameter, 
which describes the geometry of the spatial section of space-time with closed, flat and open universes 
corresponding to k = —1, 0, 1 , respectively. The coordinates r, 9 and (ft in the metric ([1]) are 'comoving' 
coordinates. The FRW models have been remarkably successful in describing the observed nature of 
universe. 

The Einstein's field equations in case of a mixture of perfect fluid and DE components, in the units 
8irG = c = 1, read as 

Rij ~^9ijR = ^ij, (2) 

where Tij = T^J 1 ' + TJ.^ is the overall energy momentum tensor with and T^ e ^ as the energy 
momentum tensors of ordinary matter and DE, respectively. These are given by 

T { f l = diag [-pM p^ J m \ p^] 

= diag [-1, u/ m ) ,J m \ w< m V m) (3) 
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and 

T {de)i = diagf-p^ p {de) J de \ P {de) ] 

= diag [-1, J de) ,J de \ u (de) ]p (de) (4) 

where p^ and p( m ^ are, respectively the energy density and pressure of the perfect fluid component 
or ordinary baryonic matter while a/ m ) = pi" 1 ) / p( m ) is its EoS parameter. Similarly, p( de ) and p( rfe ) 
are, respectively the energy density and pressure of the DE component while 

w («fc) = p (de)/ p (de) is the 

corresponding EoS parameter. 

In a comoving coordinate system, the field equations ([2]), for the FRW space-time ([T|), in case of 
([3]) and (pt|), read as 

2 ^ + ^- + \ = _ w M p M _ w (*0p(*0, (5) 
a a or 

3^ + 34 = P (m) +P (de) - (6) 
The energy conservation equation T^' ^ - = yields 

p(m) + 3^ + u (m)) p (m) H + ^(«te) + 3^ + u W) p (*) fl = Q> (7) 

where H = a/a is the Hubble parameter. 

3 Solution of Field Equations 

The field equations ([5]) and ([6]) involve five unknown variables, viz., a, u^ m \ uj^ de \ p^ and p( rfe ). 
Therefore, to find a deterministic solution of the equations, we need three suitable assumptions con- 
necting the unknown variables. 

First, we assume that the perfect fluid and DE components interact minimally. Therefore, the 
energy momentum tensors of the two sources may be conserved separately. 

The energy conservation equation ^ . = 0, of the perfect fluid leads to 

^ + 3(1+^)^)^ = 0, (8) 
whereas the energy conservation equation T^ de ^ ^ • = 0, of the DE component yields 

p(de) +3 (1+ u (de)^ p (de) H = q (g) 

Next, we assume that the EoS parameter of the perfect fluid to be a constant, that is, 
" (m) = ^ = const., (10) 

while u/^) has been allowed to be a function of time since the current cosmological data from SN la, 
CMBR and large scale structures mildly favor dynamically evolving DE crossing the PDL as discussed 
in Section 1. 

Now integration of (JH]) leads to 

/ ,N =Coa -3(iW»)) i (n) 
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where cq is a positive constant of integration. 

Finally, we constrain the system of equations with a law of variation for the Hubble parameter 
proposed by Berman [21 j . which yields a constant value of DP. The law reads as 



H = Da 



(12) 



where D > and n > are constants. In the following subsections, we discuss the DE cosmology for 
and n = by using the law (|12p . 



3.1 DE Cosmology for n^O 

In this case, integration of (|12p leads to 

a(t) = (nDt + ci)«, 
where cx is a constant of integration. Therefore, the model ([!]) becomes 

cir 2 



ds z 



-dt 2 + {nDt + ci)- 



+ r 2 (d8 2 + sm z 8 d(j)' A ) 



1 — kr 2 



(13) 



(14) 



The Hubble parameter (H), energy density (p^) of perfect fluid, DE density (p^ de ^) and EoS 
parameter (a/*)) of DE, for the model (|14p are found to be 



H = D(nDt + d) 



-i 



(15) 



p^ =c (nDt + Cl ) 



-3(l+^( m )) 



(16) 



fj \ o o -2 -3(l+u( m )) 

p (rfe) _ 3£)2^ n ^ + Ci )-2 + sk(nDt + Cl )— -c {nDt + Cl ) n . 



(17) 



(<fe) 



Jde) 



2n - 3)D 2 (nDt + Cl y 2 - k(nDt + ci)~ - c uj im \nDt + c x 



-3(1-1 



The above solutions satisfy the equation ([9]) identically, as expected. 
The spatial volume (V) and expansion scalar (6) of the model read as 

V = a 3 = (nDt + c^n, 



(18) 



(19) 



6 = 3H = 3(nDt + a)' 1 . (20) 
The density parameter r^ m ) of perfect fluid and the density parameter Q^ de ^ of DE are given by 

2n-3(l+^ m >) 



n 



(m) = = 

3# 2 3D 2 



(nDt + a) 



3F 2 L> 2 
The value of DP (g) is found to be 



3D 2 



(nDt + a 



2n-3(lW m >) 



aa 

q = g = n_ !> 



(21) 
(22) 

(23) 
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which is a constant. The sign of q indicates whether the model inflates or not. A positive sign of q, i.e., 
n > 1 corresponds to the standard decelerating model whereas the negative sign of q, i.e., < n < 1 
indicates inflation. The expansion of the universe at a constant rate corresponds to n = 1, that is, 
q = 0. Also, recent observations of SN la [1-5, 28-30] reveal that the present universe is accelerating 
and value of DP lies somewhere in the range — 1 < q < 0. It follows that in the derived model, one can 
choose the values of DP consistent with the observations. 

We observe that at t = —ci/nD, the spatial volume vanishes while all the parameters diverge. 
Therefore, the model has a big bang singularity at t = —c\/nD, which can be shifted to t = by 
choosing c\ = 0. The cosmological evolution in FRW space-time is expansionary, since the scale factor 
monotonically increases with the cosmic time. So, the universe starts expanding with a big bang 
singularity in the derived model. The parameters H, p( m \ p( de ^ and 8 start off with extremely large 
values, and continue to decrease with the expansion of the universe. The spatial volume grows with 
the cosmic time. 
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Figure 1: w (de) versus t with n = 0.05, D = 2, cq = ci = 1, u/ m) = 0. 
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The EoS parameter u/ de ) of DE asymptotically approaches to ^ — 1 provided n < 1. Thus, the 
dynamics of oj^ de ^ depends on n for sufficiently large values of t. Fig.l and Fig. 2 show the variation 
of oj( de ) during the cosmic evolution of closed (k = —1), flat [k = 0) and open [k = 1) universes for 
D = 2 and D = 6, respectively. From Fig.l, we observe that for closed and flat universes, u;( de ) starts 
in phantom region (u;( de ) < —1), crosses the PDL (u/ de ) = — 1) and finally varies in the quintessence 
region (u;( de ) > —1). However, it evolves within the quintessence region only for open universe. For 
D = 6 (see Fig. 2), u)( de ^ varies within the quintessence region throughout the evolution of the three 
universes. Thus, the nature of DE depends on the constants involved in the expression (|18p of uj^ de \ 
Also, observations predict that —1.33 < u)^ de "> < —0.79 (see, Tegmark et al. [20j). Therefore, Fig. 1 
and Fig. 2 suggest that the derived model is consistent with the observations. 

Adding (|2ip and (|22p . we get the overall density parameter 

n = ^M + ^(de) =1 + ^ nDt + c ^=^r^_ (24) 

This shows that O < 1, O = 1 and Q, > 1 according as k = —1 (closed universe), k = (flat 
universe) and k = 1 (open universe) respectively, as expected. We find that Q « 1 in the cases k = — 1 
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Figure 2: w (de) versus t with n = 0.05, D = 6, cq = ci = 1, u/ m) = 0. 



and k = 1 for sufficiently large values of t provided n < f. Thus, the model predicts a flat universe 
for large times. Since the present-day universe is very close to flat, so the derived model is consistent 
with the observations. 



1.2 
1 

0.8 
0.6 
0.4 
0.2 



0. 



n (de) 

n 



o 



0.2 



0.4 



0.6 



0.8 



Figure 3: Density parameters versus t with n = 0.05, D = 2, cq = 12, c\ = 1, a/ m ) = 0. 



Fig. 3 demonstrates the behavior of density parameters in the evolution of a flat universe (k = 0). 
We observe that initially the ordinary matter density dominates the universe. But later on, the DE 
dominates the evolution, which is probably the possible cause of acceleration of the present universe. 
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3.2 DE Cosmology for n = 

In this case, integration of (|12p yields 

a(t) = c 2 e Dt , 

where C2 is a positive constant of integration. Therefore, the model ([T]) becomes 

dr 2 



ds 2 = _ dt 2 + c 2 e 2Dt 



1 — kr 2 



+ r 2 {d9 2 + sin 2 9 dtp 2 ) 



(25) 



(26) 



The Hubble parameter, energy density of perfect fluid, DE density and EoS parameter of DE, for 
the model (|26j) are obtained as 

(27) 



H = D, 



Jm) _ -3(lW»)) -3I>(lW m >)t 



p (de) = 3D 2 + 3fcc -2 e -2^ _ CoC -3(l+u ) (-)) e -3D( 1+ .W)t ) 



CO 



(de) 



(de) 



P 



3D 2 - kc^ 2 e- 2Dt - c c- 3{1+w{m)) J^e- 3D ^ im ^ 



The above solutions satisfy the equation ([9]) identically, as expected. 
The spatial volume and expansion scalar of the model read as 



V = 4e 3Dt , 



(28) 
(29) 
(30) 



(31) 



e = 3D. 

The density parameters of perfect fluid and DE are given by 



Q(m) _ °0 C 2 



-3(l+w< m )) 



-3D(l+^( m ')t 



3D 2 



(32) 



(33) 



n (de) 



1 + 



3kc, 



3D 2 

The DP is given by 
Q=-l. 



2 e -2Dt 



c c 2 



3(l+oj( m )) 



-3D(l+w< m ))i 



3D 2 



(34) 



(35) 



Recent observations of SN la [1-5, 28-30] suggest that the universe is accelerating in its present 
state of evolution. It is believed that the way universe is accelerating presently; it will expand at the 
fastest possible rate in future and forever. For n = 0, we get q = —1 ; incidentally this value of DP 
leads to dH/dt = 0, which implies the greatest value of Hubble's parameter and the fastest rate of 
expansion of the universe. Therefore, the derived model can be utilized to describe the dynamics of 
the late time evolution of the actual universe. So, in what follows, we emphasize upon the late time 
behavior of the derived model. 
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Figure 4: versus t with D = 2, cq = c% = 1, a/ m ) = 0. 
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Figure 5: Density parameters versus t with D = 1, cq = 3, C2 = 1, Li/ m ) = 0. 



Fig. 4 depicts the behavior of the EoS parameter u/ rfe ) of DE in the closed, flat and open universes. 
We observe that for sufficiently large time, uj( de ^ ~ — 1 in each case. Therefore, the so called cosmo- 
logical constant is a suitable candidate to represent the behavior of DE in the derived model at late 
times. 

Further, at late times, we have 
p(' m ) « 

and 

« 3D 2 . 
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This shows that the ordinary matter density becomes negligible whereas the accelerated expansion 
of the universe continues with non-zero and constant DE density at late times, as predicted by the 
observations. Fig. 5 illustrates the behavior of the density parameters during the evolution of the 
universe in the derived model. It is observed that the DE component dominates the universe at late 
times and n = + fl^ » 1. 

4 Concluding Remarks 

The special law of variation for Hubble's parameter proposed by Berman in FRW space-time yields 
constant value of DP given by q = n — 1, which provides accelerating models of the universe for 
n < 1 and decelerating ones for n > 1. In the present work, we have emphasized that decelerating 
models can be described by the usual perfect fluid, while dynamics of accelerating universe can be 
described by considering some exotic type of matter such as the DE. We have investigated the role 
of DE with variable EoS parameter in the evolution of universe within the framework of a spatially 
homogeneous FRW space-time by taking into account the special law of variation of Hubble parameter. 
DE cosmologies have been discussed in two different cases, viz., power-law (n ^ 0) and exponential-law 
(n = 0). The analysis of the models reveals that the present-day universe is dominated by DE, which 
can successfully describe the accelerating nature of the universe consistent with the observations. 
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